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Abstract 

The influence of vacancies and interstitial atoms on magnetism in Pu has 
been considered in frames of the Density Functional Theory (DFT). The re- 
laxation of crystal structure arising due to different types of defects was cal- 
culated using the molecular dynamic method with modified embedded atom 
model (MEAM). The LDA+U+SO (Local Density Approximation with explicit 
inclusion of Coulomb and spin-orbital interactions) method in matrix invariant 
form was applied to describe correlation effects in Pu with these types of de- 
fects. The calculations show that both vacancies and interstitials give rise to 
local moments in /-shell of Pu in good agreement with experimental data for 
annealed Pu. Magnetism appears due to destroying of delicate balance between 
spin-orbital and exchange interactions. 

Keywords: actinide alloys and compounds, electronic band structure, 

computer simulations 
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1. Introduction 

Band structure calculations of <5-Pu predict the static magnetic order of /- 
electrons with the values of full magnetic moment 0.25-5 \xb with substantial 
impact of spin moment [l|, 0, S 5|. These results contradict to experimen- 
tal measurements of magnetic properties of non-aged Pu without impurities. 
These data indicate the absence of any ordered or disordered, static or dynamic 
magnetic moments in Pu at low temperatures 0] . 

Resent development of calculation methods allows to describe correctly the 
ground state of pure Pu in S- and model a-phase @, 0] • It has been shown in 
Ref. 0] that the delicate balance between spin-orbit (SO) and exchange interac- 
tions determines the nonmagnetic ground state in pure Pu. These interactions 
have the magnitude close to each other in actinides and its compounds and the 
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balance could be easily broken by crystal field of legands. Also, P. Soderlind [9( 
confirs the important role of SO and orbital polarization in formation of non- 
magnetic ground state of plutonium in frames of model DFT calculation. Im- 
purities like Al and Ga that are used to stabilize the fcc-phase of Pu act in the 
same way. Several groups report the presence of ordered magnetic moment in 
aged Pu-Al and Pu-Ga alloys [l(J [ll|, IT2II. The magnitude of moments is small, 
<10~ 3 p b (Ref. [n|) - 0.15 \ib (Ref. 12|), and these moments could arise 



due to distortion of crystal structure near interstitial Pu atoms and vacancies. 
Substantial drawback of the Local (Spin) Density Approximation (L(S)DA) 



is the underestimation of orbital moment [13|, [lj]. This feature leads to fail 
of the DFT in description of 4/- and 5/-mctals, since orbital moment in them 
can overcome spin one. Taking into account the Coulomb repulsion U and SO 
interactions in full matrix rotation invariant form in the LDA+U+SO method 
could improve the results. Achievement of this method is that one shouldn't 
set the exact magnetic order at the start of iterations. Both magnitude and 
direction of magnetic moment are calculated for each atom. Magnetic order and 
direction of "easy axis" are the result of self-consistent interaction procedure. 

In the LDA+U method [TBI the energy functional Elda+u depends, in ad- 
dition to the charge density p(r), on the occupation matrix n^ m , for particular 
orbital for which correlation effects are taken into account (in our case it is 
5/ plutonium orbitals). The LDA+U method in general nondiagonal in spin 
variables form was defined in Ref. [l6| 

E LDA+U [p(r), {n}} = E LDA [p(r)\ + Eu[{n}} - E dc [{n}} (1) 

where p(r) is the charge density, Elda[p(^)] is the standard LDA functional. 
The occupation matrix is defined as 

c E F 



1 r&F 

- J ImGt n >{E)dE (2) 



where G ss , (E) = (ms | (E — Hlda+u)" 1 I m's') arc the elements of the Green 
function matrix in local orbital basis set (m - magnetic quantum number, and s 
- spin index for correlated orbital) . In the present work this basis set was formed 
of LMT-orbitals from the tight binding LMTO method based on the atomic 
sphere approximation (TB-LMTO-ASA). [l7| In Eq. (TJ Coulomb interaction 
energy E'c/fjn}] term is a function of occupation matrix nf^ m , 

Eu[{n}\ = §E {mW {<m,m" | V ee \ m 1 \m'>^<X„'" 

(3) 

-(m,m" I V ee I m"',m')n 8 ^ m ,n 3 m 1, m ,„} 

where V ee is the screened Coulomb interaction between the correlated electrons. 
Finally, the last term in Eq. (TT]) correcting for double counting is a function of 
the total number of electrons in the spirit of LDA which is a functional of total 
charge density 

E dc [{n}] = ±UN(N-l)-±J H N(N-2) (4) 

where N = Tr(n^ m ,) is a total number of electrons in the particular shell. U 
and Jh are screened Coulomb and Hund exchange parameters which could be 
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determined in the constrain LDA calculations, ljj 18| The screened Coulomb 
interaction matrix elements (m. m" \ V ee \ m',m"') could be expressed via 
parameters U and Jh (see Ref. |15|). 

The functional Eq. |T]) defines the effective single-particle Hamiltonian with 
an orbital dependent potential added to the usual LDA potential 

Hlda+u = H LDA + \ ms ) V mL'(m's'\, (5) 

ms,m' s' 

V^ m , = S ss < E m » im »'{(rn,m" \ V ee \ m>, m"')n m s ,; m l + 
{(m,m" | V ee I m',m"') - (m,m" | V ee \ m'" , m'))n%„ m ,„}- 

(6) 

-U{N -\) + \j H {N -1). 

In this paper we used method LDA+U+SO which comprises non-diagonal 
in spin variables LDA+U Hamiltonian Eq. ^ with spin-orbit (SO) coupling 
term ^ ____ ^ 

Hlda+u+so = Hlda+u + Hso, 

Hso = A • L • S 

where A is a parameter of spin-orbit coupling. In LS" basis SO coupling matrix 
has diagonal (ffso)m' m as well as off-diagonal in spin variables (ffso)m< m an d 
(i?So)^I m non-zero matrix elements (in complex spherical harmonics) [? ] 

(Hso)H m = §y/(l + rn)(l-m+l)(6 m ,, m -i), 

(Hsotlm = lV( l + m)(l - m + l)(<W-i, m ), (8) 

(Hso) s m S >,m = Ams5 TO /, m 

where Im - orbital quantum numbers, spin index s = +1/2, -1/2. The pecu- 
larities of LDA+J7+SO method and its implemcnataion to the problem of pure 
Pu an several plutonium compounds were described in detales in Ref. Q . 

In the present work four different fcc-Pu supercells were investigated. Namely: 
one interstitial (IS) Pu atom in 32-atoms supercell, vacancy in 8-atoms super- 
sell, and two 32-atoms supercells with both IS and vacancy at minimal and large 
distances. Due to the presence of defects the perfect fee structure was to be 
distorted and therefore the relaxation of crystall structure for all supercells un- 
der investigation should to be taken into account. Since LMTO method doesn't 
make possible to perform structure relaxation correctly we used Classical Molec- 
ular Dynamics (CMP) with the Modified Embedded Atom Model (MEAM) by 
Baskes [20|, [2l|, [22| as interatomic potential. The MEAM is the many-body po- 
tential, i.e. interaction between a pair of atoms depends on the local structure 
(on positions of their common neighbors). The parametrization of the MEAM 
for pure plutonium and plutonium-gallium alloys was given in Ref. [2(j and 
presently the potential is widely used in CMD simulations of plutonium prop- 
erties and processes in Pu caused by self-irradiation (Refs. [2(], 21, 22, 23|, 24|). 
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Table 1: Magnetic properties calculated for 32-atoms supercell and interstitial (IS) Pu atom. 
First column - the labels of nonequivalent Pu atoms. Second column - distance between IS 
and Pu ion (A). Next four columns: number of equivalent Pu atoms in subclasses (n a toms), 
calculated values for spin (5), orbital (L), and total (J) moments. Last four columns contain 
partial contributions of / 6 configurations and jj-type of coupling for 5/ shell of Pu ion, 
effective magnetic moment (see text for explanations) and total number of /-electons. 





D. A 


H-atoms 


S 


L 


J 


f\% 


33% 


Me// 


n / 


IS 




1 


.028 


.057 


.03 


98.9 


99.2 


.146 


6.06 


Pul 


2.79 


4 


.260 


.341 


.08 


96.7 


91.7 


.234 


5.69 






2 


.065 


.059 


.01 


99.8 


97.8 


.061 




Pu2 


4.03 


4 


.216 


.310 


.09 


96.3 


93.2 


.256 


5.79 






4 


.197 


.277 


.08 


96.8 


93.8 


.238 




Pu3 


5.22 


4 


.506 


.633 


.13 


94.9 


83.5 


.271 


5.77 






8 


.380 


.472 


.09 


96.3 


87.7 


.236 




Pu4 


6.96 


4 


.750 


1.111 


.36 


85.6 


75.7 


.456 


5.73 






2 


.471 


.574 


.10 


95.9 


84.6 


.246 





Adding of IS or vacancy into initial supcreccll made Pu atoms incquivalcnt. 
That is why the different types of atoms in Tables below have additional num- 
bers (e.g., Pul etc). The relaxation of crystal structure lowers the symmetry 
again, and the new Pu classes have been divided on sub-classes (see Tab.[T|). All 
calculations of the electronic structure and magnetic properties were made us- 
ing the Tight-Binding Linear Muffin-Tin Orbitals method with Atomic Sphere 
Approximation (TB-LMTO-ASA computation scheme). In the LDA+U cal- 
culation scheme the values of direct Coulomb (U) and Hund's exchange (Jh) 
parameters should be determined as the first step of calculation procedure. It 
can be done in ab initio way via the constrained LDA calculations (l8l . Il9j |. 
In our calculations the Hund exchange parameter Jh was found to be Jh = 
0.48 eV. The value of Coulomb parameter U was set to 2.5 eV since this value 
provides the correct equilibrium volume of <5-Pu (see Ref. Q for the details). 



2. Interstitial Plutonium atom in 32-atoms supercell 

At first 32-atoms supercell of fcc-Pu with one additional Pu atom has been 
considered. The supercell has three coordination spheres around defect that is 
sufficient to describe relaxation of position of neighbor atoms. Classical molec- 
ular dynamic method was applied of describe distortion of crystal structure. 
New positions of Pu atoms in supercell were used in further calculation of elec- 
tronic structure. Adding of one additional Pu atom lowers the symmetry of 
the cell. Four new inequivalent classes of plutonium belonging to four differ- 
ent coordination spheres around IS arouse. Moreover Pu atoms within each 
new class become inequivalent due to different local neighborhood. To take 
in to account this lowering of symmetry no symmetrisation was applied in our 
electronic structure calculation. Since no additional symmetry conditions was 
imposed on electronic subsystem magnitude of local moments on Pu sites and 
their directions can be arbitrary and correspond to the minimum of total energy. 

The LDA+U+SO calculations for metallic Pu in S phase gave a nonmagnetic 
ground state with zero values of spin S, orbital L, and total J moments [1, @|. 
Our calculation made for 32-atoms supercell with one IS shows that small local 
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magnetic moments develop on Pu sites. The magnitude of local moment depends 
on the distance between center of distortion (IS) and corresponding Pu site. 
We argue that this appears because of the break of balance between SO and 
exchange interactions due to the relaxation of crystal structure. The results are 
presented in Tab. [TJ The partial contributions of / configuration and jj-type 
of coupling to the final state could be calculated in the following way. Total 
moment value is the same in both coupling schemes (jj or LS): J = for / 6 
and J = 5/2 for / 5 . If there is a mixed state (1 — x) ■ f 6 + x ■ f 5 then x can be 
defined as x = J/2.5. Spin S and orbital L moment values for f e configuration 
are equal to zero in jj coupling scheme and S = 3, L = 3 in LS coupling 
scheme. For f b configuration they arc £ = 5/14 « 0.36, L — 20/7 ~ 2.86 in jj 
coupling scheme and S* = 5/2,L = 5in LS coupling scheme. One can define 
a mixed coupling scheme with a contribution of jj coupling equal to y and of 
LS coupling to (1 — y), correspondingly. In final state the calculated values of 
orbital and spin moments will be: 

L = x ■ (2.86 ■ y + 5 • (1 - y)) + (1 - x) ■ (0 • y + 3 • (1 - y)), (9) 

S = x ■ (0.36 • y + 2.5 • (1 - y)) + (1 - x) ■ (0 • y + 3 • (1 - y). (10) 

These formulas allow to determine the value of coefficient y. An effective 
paramagnetic moment obtained from susceptibility measurements using Curie- 
Weiss law can be calculated as 

fleff = 9 ■ \/ J ■ {J + 1) • Ms- (11) 

The problem is to define Lande g-factor which can be calculated for pure f 5 
and f e configurations in LS or jj coupling schemes. As for / 6 configuration 
total moment J = 0, one needs to calculate g- factor for / 5 configuration only. 
For ground state of / 5 configuration in jj coupling scheme Lande factor is gjj 
= 6/7 ~ 0.86. In LS coupling scheme its value is qls = 2/7 « 0.29. As the 
latter value is nearly three times larger than the former, gjj and gLS can give 
only upper and lower limits of g-factor for the case of intermediate coupling. 
We could calculate weighted value of effective moment using relative weights of 
LS- and jj-couplings obtained from Eqs. [101 [HI 

The magnetic order of Pu ions is set arbitrary at the beginning of itera- 
tion process. Final directions of local moments are calculated in accordance 
with minimum of total energy at the end of self-consistence loop. Since the 
long range order is set up in used calculation method, some fcrrimagnctic order 
arouse as an artifact. This order resembles the antifcrromagnetic (AFM) one of 
A-type. In Fig. [T] resulting directions of total moment are shown with red and 
green colors. Pu atoms positioned in the first coordination sphere to IS as well 
as IS itself have the smallest magnetic moments. They don't differ significantly 
from pure <5-Pu that is nonmagnetic. The values of local moments grow up with 
increasing distance between IS and corresponding site. The largest total mo- 
ment developed for Pu4 ion positioned in the center of supercell (large red sphere 
in Fig.[T]). This ion has the largest distance to IS. The average value of the effec- 
tive moment in the case of IS Pu atom in 32-atoms supercell is Heff ~0.26 (1b- 
Numbers of /-electrons (see Tab. [TJ last column) differ from those calculated for 
not distorted fcc-Pu (which has 5.74 /-electrons) but not significantly except 
the case of interstitial atom. The later has the largest occupation number in all 
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Figure 1: (Color online). Non-collinear order obtained for 32-atoms supercell with IS. Purple 
spheres denote IS. Green and red spheres are Pu atoms with opposite directed total moments. 
The radii of spheres are proportional to the magnitude of correspondent magnetic moments. 
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Table 2: Magnetic properties of Pu ions calculated for 8-atoms supercell with one vacancy. 
Second column - distance between vacancy and Pu ion (A). See also the caption of Tab. ^ 





D. A 


^atoms 


S 


L 


J 


f\% 


33% 


Me// 


n / 






2 


.651 


.817 


.166 


93.4 


78.7 


.298 




Pul 


3.27 


2 


.518 


.656 


.138 


94.5 


83.1 


.283 


5.71 






2 


.515 


.652 


.137 


94.5 


83.3 


.283 




Pu2 


4.63 


1 


.356 


.449 


.094 


96.3 


88.5 


.243 


5.70 



considered structures (see tables below) . The large number of /-electrons (close 
to 6) obtained in the present calculation disagrees wit h p revious experimental 
and theoretical estimations that give 5.1-5.2 electrons [25j. Such difference be- 
tween theoretical results originates from dissimiliar band structure calculation 
methods. Since TB-LMTO-ASA scheme use artifically large overlaping atomic 
spheres these numbers should be considered only to compare differnt clases of 
plutonium atoms with each other. We have verified our results and run sev- 
eral calculations with diffcrnet radii of Pu atoms filling of empty space in the 
primitive cell with empty spheres (pseudo-atom without core states) . Distorted 
supercell always becomes magnetic with the same order. The magnitude of local 
magnetic moments depends slightly on atomic radius. It rises up with increas- 
ing of the later. For the sake of simplicity we have choosen the same radii 3.41 
a.u. for all Pu atoms to be able to compare their magnetic moments. Artifical 
overlapping of atomic spheres in all considered supercells never overcomes 13% 
which is critical TB-LMTO-ASA value. 

3. Vacancy in 8-atoms supercell 

Another type of defects appearing in fresh Pu during first several years [26| is 
vacancy. Small supercell consisting of 8 Pu atoms was considered. One Pu atom 
was removed from its position in the supercell and after relaxation of crystal 
structure this empty space was artificially filled with empty sphere. Unfortu- 
nately, 8-atoms supercell is not sufficient to describe correctly the relaxation 
of crystal structure within the MEAM. Shifts of Pu atoms were obtained to 
be negligible. Nevertheless removing of one atom from the supercell lowers the 
symmetry of crystal, since the local neighborhood of plutonium atoms becomes 
diffcrnet. Note, that this turns out enough to give rise to the local moments on 
Pu sites. Opposite to IS action, the vacancy affects much larger on Pul that 
form the first coordination sphere of it (Tab. [2]). Magnetic moments of Pu2 
atoms that belong to the 2-nd coordination sphere are smaller. 

The average value of the effective magnetic moment in the supercell with 
vacancy is ~0.28 fiB- In contrast to the case of IS in 32-atoms supercell, the 
resulted magnetic order is the analogue of C-type AFM (see Fig. 0). 

These results prove that both types of defects call local magnetic moment 
on Pu atoms due to distortion of fee structure or even lowering of symmetry. 
Different types of defects affect magnetism in Pu in different ways: IS calls the 
largest magnetinc moment on an atom at large distance whereas vacancy affects 
mostly its nearest neighbors. Different types of defects result also in different 
types of AFM order. Simultaneous effect of IS and vacancy could also give rise 
to the local moments and lead to more complicated pattern of Pu ions magnetic 
order. 
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Figure 2: (Color online). Non-collinear order obtained for 8-atoms supercell with vacancy. 
Vacancy is shown as black spheres. Green and red spheres are Pu atoms with opposite directed 
total moments. The radii of spheres are proportional to the magnitude of correspondent 
magnetic moment. 
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Table 3: Magnetic properties of Pu ions calculated for 32-atoms supercell with one vacancy 
and IS at minimal distance. Sec also caption of Tab. IT1 





n a 


Hatom 




T 


T 
J 


J ,/o 




Me// 


n/ 


Id 




1 
1 


U.UOo 


U.Uou 


n D9Q 
U.uzy 


Q8 Q 

yo.y 


QR 9 

yo.z 


U. 140 


p. i 

0. 1 






o 
Z 


u.zou 


U.o4o 


u.uy i 


yo. i 


09 n 
yz.u 


n 9^q 
u.zoy 




Pul 


4.02 


9 
z 


n 1 Q 1 ^ 
u. lyo 


n 9fi7 

U.Zu t 




Q7 1 
y / . i 


yo.o 


fl 99^ 
u. zzo 


5.61 


rUZ 


9 70 
Z. / U 


i 
1 


u.uuy 


U.UIO 


n nnfi 
u.uuu 


QQ 8 

yy .o 


QQ 7 

yy. / 




f! 99 
u.zz 






2 


0.485 


0.583 


0.098 


96.1 


84.1 


0.237 




Pu3 


5.10 


2 


539 


645 




Q5 5 


89 6 


0.252 


5.71 






2 


0.230 


0.313 


0.083 


96.7 


92.7 


0.239 




Pu4 


6.91 


2 


0.021 


0.038 


0.017 


99.3 


99.4 


0.111 


5.75 






2 


0.165 


0.212 


0.047 


98.1 


94.7 


0.180 




Pu5 


4.09 


2 


0.122 


0.155 


0.033 


98.7 


96.1 


0.153 


5.71 






2 


0.057 


0.081 


0.024 


99.1 


98.2 


0.132 




Pu6 


2.85 


2 


0.088 


0.117 


0.029 


98.8 


97.2 


0.145 


5.78 


Pu7 


5.18 


2 


0.123 


0.156 


0.034 


98.7 


96.1 


0.154 


5.75 


Pu8 


6.91 


1 


0.069 


0.077 


0.008 


99.7 


97.7 


0.074 


5.75 


Pu9 


5.35 


2 


0.220 


0.277 


0.057 


97.7 


92.9 


0.196 


5.77 


PulO 


6.99 


2 


0.343 


0.448 


0.105 


95.8 


88.9 


0.261 


5.78 






2 


0.069 


0.095 


0.027 


98.9 


97.8 


0.140 




Pull 


5.24 


2 


0.150 


0.199 


0.049 


98.0 


95.2 


0.186 


5.79 



4. Vacancy and Interstitial Plutonium at minimal distance 

Since IS and vacancy affects on the magnetism in Pu in different ways, we 
can expect that their simultaneous influence could also give rise of local moment, 
and produce some complicated magnetic pattern. 

As the first step, 32-atoms supercell with IS and vacancy at minimal distance 
has been investigated. Relaxation of the supercell was made with use of molec- 
ular dynamics within MEAM. Both removing one Pu atom from its site and 
relaxation lowers the symmetry of considered supercell and one calss of Pu has 
been devided in 11 classes. Six of them have 2 subclasses (see Tab. [3]). Small 
local moments also developed at Pu sites of relaxed supercell. Magnitude of 
moments and other results of LDA+U+SO calculation are presented in Tab. [3] 
Two different values of moments for one type of Pu atoms is the consequence 
of lowering of the symmetry of supercell due to orbital polarization. Mutual 
action of IS and vacancy decreases dispersion of magnitudes on differnet sites. 
The averaged value of effecive moment at Pu atom is 0.18 ^b- 

Three types of Pu atoms, PulO, Pul, and Pu3, have the largest magni- 
tudes of magnetic moment, 0.261 fis, 0.259 fiB, and 0.252 (is, correspondingly 
Atoms of type Pul and Pu3 are the nearest to the vacancy (except IS) and have 
hence the largest magnetic moments in agreement with results of our previous 
calculation for vacancy in 8-atoms supercell (Sec.[3]). Pu6 atoms belong to the 
first coordination sphere of vacancy but posses much smaller moment, ~.14 /xb. 
This atom is positioned in the first coordination sphere of IS and, in agreement 
with our results for one IS in 32-atoms supercell, IS suppresses magnetism on 
Pu6 atom. Finally, PulO has a sizeable value of magnetic moment but it is 
smaller then that at Pu atom in the center of 32-atom supercell with single IS. 
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This could be explained by the action of vacancy that calls large local moments 
near itself and suppresses the magnetism on far standing atoms. 

Simultaneous effect of IS and vacancy results in more complicated canted 
AFM pattern, which could not be identified with any standart type. Calculated 
canted AFM order is presented in Fig. [3l 

5. Vacancy and Interstitial Plutonium at large distance 

Finally, we made the same calculation for the 32-atom supercell contained 
IS and vacancy at the large distance. As well as in previous cases relaxation 
of crystal structure has been made within MEAM before band structure cal- 
culation. Considered defects also lower the symmetry and 17 new Pu classes 
arouse. The values of moments and contribution of coupling types and electronic 
configuration are presented in Tab. 2] 

Like our calculation for IS and vacancy at minimal distance, local moments 
developed at all Pu atoms. The magnitude of moments depends on the dis- 
tances to both IS and vacancy. Incommensurate magnetic order with strong 
noncollinearity was obtained for this type of defects position (see Fig. [4]) . The 



b 




Figure 3: Ferrimagnetic order obtained for 32-atoms supercell with IS and vacancy at minimal 
distance. Black spheres denote vacancy, purple spheres - IS. Green and red spheres are Pu 
atoms with opposite signs of z-component of magnetic moment. The length of arrows is 
proportional to the magnitude of correspondent effective magnetic moment. 
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Table 4: Magnetic properties calculated for 32-atoms supercell with one vacancy and IS at 



large distance. See also caption of Tab. ^ 





— 77 5 

D,A 


diatoms 


S 


L 


J 


/ 6 ,% 


33,% 


Me// 


n / 


IS 




l 


.052 


.083 


.032 


98.7 


98.4 


0.153 


6.12 


Pul 


4.26 


4 


.144 


.176 


.032 


98.7 


95.4 


0.149 


5.70 


Pu2 


2.79 


2 


.103 


.141 


.038 


98.5 


96.7 


0.164 


5.78 


Pu3 


5.01 


1 


.235 


.333 


.098 


96.0 


92.6 


0.261 


5.73 


Pu4 


6.67 


2 


.152 


.188 


.035 


98.6 


95.1 


0.157 


5.72 


Pu5 


3.77 


4 


.268 


.361 


.093 


96.3 


91.4 


0.250 


5.70 


Pu6 


3.11 


1 


.065 


.091 


.025 


99.0 


97.9 


0.136 


5.73 


Pu7 


5.25 


2 


.031 


.039 


.007 


99.7 


99.0 


0.072 


5.74 


Pu8 


7.12 


1 


.173 


.219 


.046 


98.2 


94.4 


0.178 


5.71 


Pu9 


2.63 


1 


.077 


.128 


.051 


98.0 


97.7 


0.194 


6.01 


PulO 


4.90 


2 


.226 


.276 


.050 


98.0 


92.7 


0.182 


5.75 


Pull 


6.83 


1 


.169 


.276 


.107 


95.7 


94.8 


0.280 


5.73 


Pul2 


2.79 


2 


.020 


.031 


.012 


99.5 


99.4 


0.093 


5.82 


Pul3 


4.87 


2 


.139 


.181 


.042 


98.3 


95.5 


0.172 


5.78 


Pul4 


7.29 


2 


.073 


.131 


.058 


97.7 


97.8 


0.208 


5.75 


Pul5 


5.23 


2 


.239 


.303 


.065 


97.4 


92.3 


0.208 


5.77 


Pul6 


5.42 


1 


.164 


.228 


.064 


97.5 


94.8 


0.212 


5.78 


Pul7 


5.02 


1 


.245 


.312 


.068 


97.3 


92.1 


0.213 


5.79 



mechanism of the formation of magnetic moments was described above. The 
averaged magnetic moment on Pu is 0.179 [Ib- 



6. Conclusions 

Band structure calculations have been run for 4 supercells contained IS, va- 
cancy, and both IS and vacancy at small and large distances. For the supercell 
with one IS ferrimagnetic order close to A-type AFM was obtained. The mag- 
inudes of local moments are 0.06-0.46 \ib and the averaged moment is 0.26 (1b- 
The largest value of magnetic moment has atoms at the largest distance from 
IS. Ferrimagnetic order close to C-type AFM was obtained for 8-atoms super- 
cell with vacancy. This type of defect calls the largest moment on Pu atoms 
in first coordination sphere. Magnetic moments obtained in 32-atoms super- 
cell with both IS and vacancy have smaller dispersion of magnitude 0.1-0.3 /is 
and smaller averaged moment ~0.18 fiB- Simultaneous action of these defects 
results in incommensurate magnetic order with strong non-collinearity. Never- 
theless the long range order obtained in the present work should be considered 
as an artifact of computation method. Our results indicate that sort range or- 
der could appear due to defects in fee Pu. And the type of such order depends 
strongly on the distance to the corresponding defect. The implementation of the 
LDA+DMFT method is necessary to describe more accurately the magnitudes 
of local moment in paramagnetic phase of Pu. 

Results of calculation explain presence of magnetic moment in aged Pu sam- 
ples and agree well with experimental data 11, 12, (|. 



Support by the Russian Foundation for Basic Research under Grant No. 
RFBR-07-02-00041 is gratefully acknowledged. 



11 




Figure 4: Ferrimagnetic order obtained for 32-atoms superccll with IS and vacancy at large 
distance. See also caption of Fig. [3] 
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